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A COMPARISON OF MEASURED AND PREDICTED TEST FLOW IN AN 
EXPANSION TUBE WITH AIR AND OXYGEN TEST GASES 


Kenneth V. Haggard and William K. Goad 
Langley Research Center 

SUMMARY 

Simultaneous time-resolved measurements of temperature, density, pitot pressure, and 
wall pressure in both air and O 2 test gases have been obtained in the Langley pilot model 
expansion tube. These tests show nonequilibrium chemical and vibrational relaxation signifi- 
cantly affect the test-flow condition. The use of an electromagnetic device to preopen the 
secondary diaphragm before the arrival of the primary shock wave resulted in an improvement 
in the agreement between the measured pitot pressure and the value inferred from measured 
density and interface velocity. Boundary-layer splitter plates used to reduce the wall boundary 
layer have shown that this disagreement in the measured and inferred pitot pressures is not, 
as speculated in NASA TN D-7273, a result of boundary-layer effects. 

INTRODUCTION 

For several years a research program has been conducted at Langley Research Center to 
determine a range of operating conditions for which the expansion tube flow is suitable for 
hypervelocity testing. Experimental studies performed in the Langley pilot model expansion 
tube are described in references 1 to 4, and several idealized theoretical flow models are 
described in references 5 to 10. The complexity of the expansion tube flow cycle is such 
that, in spite of this extensive long-term effort, it has not been possible to develop a flow 
model that can consistently predict the test-flow conditions or explain the relationships between 
the measured test-flow quantities. 

The study described in reference 1 used simultaneous measurements of wall pressure, 
pitot pressure, density, and the vibrational temperature of to demonstrate that chemical 
and vibrational nonequilibrium processes were probably important in the flow cycle. These 
results are in qualitative agreement with the predictions of references 8, 9, and 10. However, 
the temperature and density measurements were dependent on absorption measurements made 
across the tube diameter and were thus influenced by the wall boundary layer. Self-consistency 
checks among the four measurements and estimates of the wall boundary-layer thickness 
obtained from reference 4 indicated there could be an error of as much as 20 percent in the 
measured density due to the boundary layer. 



The present study was undertaken to determine the effects of nonequilibrium chemical 
and vibrational relaxation on an expansion tube test flow. These tests were conducted with 
both pure O 2 and air as a test gas, and the range of test conditions was extended beyond 
that of reference 1 to better determine the appropriate flow model for the facility. These 
studies used a set of splitter plates in the test section to reduce the effects of the wall 
boundary layer on the density and temperature measurements. 

Some experiments used an electromagnetic device to open the secondary diaphragm prior 
to the arrival of the shock wave. This device could conceivably improve the quality of the 
test flow by reducing the contamination from the diaphragm material and by eliminating the 
need for the test flow to supply the momentum necessary to open the diaphragm. 

SYMBOLS 


absorption parameter 
a sound speed, m/s 

F mass fraction of O 2 

H enthalpy, m^/s^ 

*X intensity at X (arbitrary units) 

*X,o unabsorbed intensity at X (arbitrary units) 

k X absorption coefficient at X, cm - * 

L distance from secondary diaphragm to test section, m 

1 length, cm 

M molecular weight of initial test gas 

p pressure, Pa 

p t pitot pressure, Pa 

p pitot pressure defined by equation (9), Pa 

t,c 


2 



wall pressure, Pa 


universal gas constant, 8.31 J/mole-K 
temperature, K 

vibrational temperature of O^, K 
time 

time after arrival of test flow at test section, s 
residence time of fluid particle behind primary shock, 
flow speed, m/s 

speed of secondary interface, m/s 
speed of primary shock wave, km/s 
distance, m 
compressibility factor 
ratio defined in equation (3) 



wavelength, nm 
density, kg/m^ 

scale depth, cm (see eq. (5)) 
density of O^, kg/m^ 


Subscripts: 


0 standard conditions (0.101325 MPa; 273.3 K) 

1 initial state in intermediate section 

2 state behind primary shock 

5 state of test-gas flow 

6 state behind reflected shock at secondary diaphragm 

10 initial state in acceleration section 

DESCRIPTION OF FACILITY 

The Langley pilot model expansion tube was a pilot facility designed to study the 
feasibility of using an expansion tube to produce a short duration hypersonic flow duplicating 
conditions encountered by an entry vehicle. This pilot model facility was terminated when 
the Langley 6-inch expansion tube at Langley Research Center became operational (see ref. 11). 
The pilot model tube was a 9.55-cm-diameter stainless steel tube consisting of three sections: 
a driver section 1.55 m long, a 10-m intermediate section, and a 9.75-m acceleration section. 
The driver and intermediate sections were separated by a 3.17-mm-thick steel diaphragm. The 
intermediate and acceleration sections were separated by either a 0.0063 5-mm Mylar diaphragm 
or a 0.4-mm-thick aluminum sheet. 

An idealized distance-time diagram of the flow cycle is shown in figure 1(a). In order 
to start the cycle, the intermediate chamber was filled with the test gas (state Q)) and the 
acceleration section was filled with helium (state (£§))• The driver section was then pressur- 
ized with hydrogen (state until the first diaphragm ruptured. At this time the primary 

shock wave was formed; as the shock wave moved down the intermediate section, it com- 
pressed and heated the test gas to state (5) . The shock was reflected when it reached the 
Mylar diaphragm, and the increased pressure caused the diaphragm to rupture. The details 
of this process are shown in figure 1(b), where it has been assumed that the diaphragm opens 
instantaneously a short time after the shock reflection. A secondary shock wave is formed 
and moves down the acceleration section when the Mylar diaphragm ruptures, and a second 
expansion fan provides the transition from the reflected shock region (state (6) ) to the final 
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test flow ^state(J)). For the conditions at which these tests were made, the test flow is 
ended by the arrival of the tail of the expansion fan at the test section. 

INSTRUMENTATION 

The instrumentation used to measure pressure and wave speeds is described in refer- 
ences 1 and 2. The speed of the shock waves was measured by using a combination of 
pressure transducers, ion gaps, and microwave reflection. The interface speed was usually 
determined by microwave reflection with a back up of photodetectors and pitot pressure mea- 
surements. A description of the microwave reflection technique is contained in reference 12. 

A sketch of the arrangement of the equipment used to measure density and temperature 
at the test section is shown in figure 2(a). The ultraviolet source has a capillary tube through 
which hydrogen was continuously pumped. An artificial transmission line was connected across 
the lamp and a spark series gap. A synchronized pulse triggered by the primary shock breaks 
down the spark gap so that the artificial transmission line discharges through the lamp, thereby 
producing an intense source of vacuum ultraviolet radiation that lasts for more than 2 ms. 

The light from the lamp passes through the test section of the tube via two CaF^ windows 
and into a 0.5-m vacuum ultraviolet spectrograph. The intensity of the diffracted beam was 
measured at 147 and 16U nm over a bandwidth of 0.8 nm by sodium salicylate coated 
1P21 photomultiplier tubes. The output of each phototube was displayed on an oscilloscope 
and photographed. 

A sketch of the test section with the splitter plates in position is shown in figure 2(b). 
The plates were 36.5 mm wide and 53 mm long, with the window 25 mm behind the leading 
edge. The plates were spaced 27 mm apart. Wall pressure was measured just upstream of 
the light beam, and the pitot pressure was measured just downstream of the beam and on the 
center line of the tube. For the conditions investigated, the splitter plates did not affect the 
pitot pressure measured by the probe during that part of the flow covered by these studies. 

It is estimated that there is typically a ±20-percent error in the temperature measure- 
ment and a ± 10-percent error in the density and pressure measurements. 

DIAPHRAGM OPENER 

The secondary diaphragm is generally opened by the increase in pressure behind the 
reflected shock in the test gas. Thus the diaphragm is subjected to the high temperatures of 
the reflected shock region; and although the exposure time is short — a few microseconds — 
surface contamination on the diaphragm is transferred to the test gas. In addition, the test 
gas loses momentum to the diaphragm during the opening process. In an attempt to reduce 
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these problems without introducing other major detrimental effects, an electromagnetic device 
was used to open the diaphragm prior to the arrival of the primary shock wave. 

A sketch of the diaphragm-opener coil is shown in figure 3. It consisted of an insulate! 
aluminum coil with an inner diameter of 9.55 cm and an outer diameter of 14.68 cm. The 
two parallel plate prongs served to connect the coil to a high voltage capacitor bank. The 
installation of the coil in the expansion tube is shown in figure 4. The prongs extended out 
of the tube through a double O-ring seal, and the body of the coil was supported inside the 
tube by a fiberglass liner. Upstream of the diaphragm location, the inner diameter of the 
liner matched that of the tube; however, downstream of the diaphragm location, the inner 
diameter opened up to the outer diameter of the coil and then tapered down to the tube 
diameter. No studies were made of the effects on the flow of the area change in the tube 
downstream of the coil. A 0.4-mm-thick aluminum sheet was used as a diaphragm. The 
center of the diaphragm was scribed with a cross to half thickness, and the diaphragm was 
clamped into position so that it made contact with the downstream side of the coil. 

When the capacitor bank was discharged through the coil, the electromagnetic field about 
the coil generated eddy currents in the diaphragm; consequently, the diaphragm was forced 
away from the coil, causing it to tear along the scribes. Once the tabs had been given the 
initial momentum, the fact that they were clamped at the outer edge caused them to swing 
outward into the enlarged area just downstream of the coil where they were out of the direct 
flow. So long as an excessive amount of energy was not transferred to the diaphragm, the 
tabs did not tear loose and the opened tabs remained in the enlarged area throughout the 
rest of the flow cycle. 

A series of photographs of the opening process is shown in figure 5. These photographs 
were taken during a bench test of the opener. The bank voltage was 5.3 kV and the energy 
was 10.8 kJ. The first photograph (fig. 5(a)) shows the initial tearing of the diaphragm along 
the scribes. At this time the diaphragm is still conducting a large eddy current, and the 
sudden tearing of the diaphragm causes the arcing across the tears between the tabs. The 
second and third photographs (figs. 5(b) and (c)) were taken at 10-ps intervals and show the 

arcing peak and start to diminish. The other photographs show the continued opening of 

the tabs. 

For a large number of bench tests and for eight expansion tube tests, the diaphragm 
opener performed well mechanically. However, on the ninth expansion tube test, the trigger 
signal failed to discharge the bank and the diaphragm was forced open by the reflected shock 

wave. The coil seemed to have survived this trouble with no damage; however, on the next 

test the insulation on the coil broke down and the coil was shorted. This breakdown 



destroyed the coil and ended the testing with the opener. The authors feel it is probable 
that the coil insulation had been damaged in the eighth run since the insulation had withstood 
repeated bench tests under voltages as high as 8 kV. 

Because of the limitations in the facility, it was not possible to install the microwave 
measurement device in the tube with the diaphragm opener in use. Thus, the interface 
velocities for these runs were obtained from phototube records and the pitot pressure mea- 
surements. Although no tabs were tom off during these tests, the possibility was such a con- 
stant worry that the splitter plates were removed from the tube during the diaphragm-opener 
runs. 


DATA REDUCTION 


The vibrational temperature and density of C >2 were determined from the absorption 
measurements by using the method described in references 13 and 14. Basic data and a 
theoretical treatment are to be found in references 15 and 16. With this technique the trans- 
mission of the flow is measured at two wavelengths (147 and 160 nm). These wavelengths 
were selected because the absorption coefficient at 147 nm is sensitive to temperature whereas 
the absorption coefficient at 160 nm is relatively insensitive to temperature. Thus, since the 
absorption is only a function of temperature, density, and wavelength, the temperature and 
density can be determined. This procedure will become clear by considering the equations 
governing the absorption process. 

The attenuation of the ultraviolet radiation was assumed to follow Beer’s law. Thus, 


where Pq 
tion path. 


*X (t) = 1 X,o (t)exp 



p(t)k^(t) 


dl 


( 1 ) 


is the absorber density in this case J and the integral is taken over the absorp- 

The absorption parameter A. is defined by the relation 


A x (t) = log e 


■x<*> 


= ^ d> 


( 2 ) 


and 


a, the ratio of Aj 47 


a(t) = 


A 147 (t) 

A 160 (t) 


to A 160 , is 

f P(t)k I47 (t) dl 

jf p ( t)k 160 (t) dl 


where the subscripts 147 and 160 are values of X in nanometers. 


(3) 
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The data-reduction model for the density and vibrational temperature of includes 
the following three assumptions: 

(1) A valid vibrational temperature exists in the flow (i.e., the vibrational energy levels 

have a Boltzmann distribution) 

(2) Absorption by species other than is negligible 

(3) Flow quantities are uniform along the absorption path. 

With these assumptions, a becomes 


k I47O) 


(4) 


Plots of kj 4 y, kjgQ, and a as a function of temperature are shown in figure 6. These 
curves have been taken from reference 14. 

Once a is obtained from measured values of A^, T y can be found from 
figure 6. The scale depth p^ is defined by 

f p(t) ^x(f) 

s J V d ’ = w <5) 

The O 2 density is then obtained from 

p o 2 ,o 

p o 2 (t) = T p d(t> 


where t is the path length across the flow. 

When each of the values p w , P-j-, Pq , T y q , and Vj have been measured, they 
can be cross-checked by using the equation of^ state and by assuming that the flow is in 
thermodynamic equilibrium. The flow density is then obtained from 

P (t) = f P G (t) (7) 
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where the mass fraction F of in the flow can be calculated from the measured pressure 
and temperature. Thus, the values of each measured quantity can be compared with the 
value obtained from the other measured quantities and from the relation 

_ - 


where Z is assumed to be unity. If, in addition, it is assumed that the flow velocity is 
given by the interface velocity, then the pitot pressure is approximately given by 

P t ,C = P V i 2 <»> 


This relation (eq. (9)) is accurate to within 3 percent for the conditions tested and is valid 
even if the flow is not in thermodynamic equilibrium. Thus, if A is defined as 


A = 



J o 


V.' 


F Pt 


( 10 ) 


and is set equal to unity, the mass fraction of in the flow can be found from 


F = 


O 


When a value of A is given, it has been obtained by assuming that 
value in the initial test gas ^state ©. fig- i). 


do 

F is defined by its 


A detailed discussion of the measurement errors in the vibrational temperature is given 
in reference 1. In general, errors of 20 percent in T v and 10 percent in density and 
pressure are expected. Measurements of the shock and interface speeds are accurate to 1 or 
2 percent. 


THEORETICAL FLOW MODELS 

Wave diagrams similar to those in figure 1 can be used to construct mathematical models 
of the expansion tube flow cycle (see, for example, refs. 5 to 10). Because of the complex 
nature of the expansion tube flow, each of these models involves major simplifying assump- 
tions. For example, they all neglect viscous effects and interface mixing. The model described 
in reference 7 treats the finite opening time of the secondary diaphragm, and the models 


St 
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described in references 8 to 10 treat finite-rate chemical and vibrational relaxation in part of 
the flow cycle. 

There are five flow models that seem to provide good upper and lower bounds for the 
measured flow quantities. These models (ISPG, ISEE, RSEE, ISFE, and RSFE) are described 
as follows and are discussed in more detail in references 1 and 17: The ISPG model assumes 

the test gas is a perfect gas. The flow quantities are calculated by using the flow diagram in 
figure 1(a) (no reflected shock); the expansion process is assumed to be isentropic. The ISEE 
model, like the ISPG model, assumes there is no reflected shock at the secondary diaphragm 
and that the expansion process is isentropic. However, the test gas is assumed to be real, 
and the equilibrium properties of the fluid are used throughout the flow cycle. The RSEE 
model assumes there is a reflected shock upstream of the secondary diaphragm and that the 
entire test flow is processed by this shock. The expansion is assumed to be isentropic, and 
equilibrium flow properties are used throughout the entire calculation. This RSEE model 
neglects any effects of the interaction of the head of the expansion fan with the reflected 
shock. The ISFE and RSFE models are similar to the ISEE and RSEE models except that 
the ISFE and RSFE models assume the flow becomes frozen upon entering the expansion fan. 
That is, the chemical composition and vibrational energy are fixed at the equilibrium values 
just upstream of the head of the expansion fan, and the chemical composition and vibrational 
energy remain constant during the rest of the flow cycle. 

These five models neglect the finite rates of chemical and vibrational relaxation; however 
the frozen-expansion and equilibrium-expansion processes should provide limits on this non- 
equilibrium problem. The only attempt to treat the nonequilibrium character of the flow 
theoretically is described in references 8, 9, and 10. The model described therein assumes the 
flow is in equilibrium until it enters the expansion fan and only then allows for finite-rate 
chemical and vibrational relaxation. This assumption breaks down for particle paths near the 
diaphragm since the flow along these paths cannot reach equilibrium before entering the 
expansion fan. The severity of this breakdown can be seen by examining the nonequilibrium 
flow behind the primary shock for a typical test condition. For such a case the mole fraction 
of oxygen atoms as a function of time is shown in figure 7. The test gas is pure and 
the shock speed is 2956 m/s, and the initial temperature is 300 K. The equilibrium mole 
fraction of O is calculated to be 0.115. The abscissa is the residence time t res of the 
fluid particle between the shock into state (T) and the head of the expansion fan (this assume 
there is no reflected shock). By assuming that states ( 2 ) and (?) in figure 1(a) are uniform, 
t res can be related to the test time t^ of that fluid particle. The test time t<- is the 
time after the arrival of the secondary interface at the test section. Thus 


Ves 


55 >1 V S,1 


p 2 a 2 t 


( 12 ) 
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For the conditions of this test the test flow is ended by the arrival of the tail of the expan- 
sion fan at the test section. Thus the maximum test time t. is 


*5, max 



(13) 


For this same test, the maximum residence time is 109 ps for the ISEE model but 23.5 /us 
for the ISFE model. These times are so marked by the circles in figure 7. Clearly, the 
flow does not reach equilibrium behind the primary shock; and if a nonequilibrium calculation 
of the flow cycle is to be undertaken, it must start with the primary shock, not with the 
expansion fan. 

All of the calculations for the various theoretical models used in this report were made 
with the use of the computer code described in reference 1 7. The nonequilibrium data for 
Figure 7 were obtained from the computer code of reference 18. 

TEST CONDITIONS 

The initial pressures, the test-gas composition, and the measured primary shock and 
secondary interface speeds are listed in table 1. In all cases the driver section contained 
hydrogen, and the acceleration section either contained helium or was evacuated. The initial 
test-gas pressure Pj varied from 2933 to 333 Pa and the acceleration section pressure 
Pjq, from 133 to less than 1 Pa. Primary shock speeds ranged from 4002 to 2718 m/s, 
and the interface velocity varied from 6458 to 4938 m/s. 

RESULTS AND DISCUSSION 
Time Variations 

Because of the broad range of conditions tested, it is not possible to pick one condition 
that is representative of the entire range covered in this report; however, there is one condition 
that produces test results for which the measured values of p w , p^, p, T v o 2 > ar *d Vj 
are in good agreement with each other when equations (8) and (9) are used to compare the 
flow quantities. This condition is run 3676 in table 1. Time histories of an example of 
the run at this condition are shown in figure 8. The test gas is at an initial pressure 

of 2933 Pa; the primary shock speed is 2955 m/s, and the flow was expanded to 6101 m/s. 
The measured values of p w , p t , p, and T v q 2 are denoted by the continuous curves and 
are plotted as a function of the test time t^. For each flow quantity, equation (8) or (9) 
has been used to infer the corresponding value based on the other measured values. These 
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inferred values have been averaged over 25-ps intervals to reduce the scatter of the experimen- 
tal data and are denoted in figure 8 by the circular symbols. 

The comparison between the measured and calculated values is quite good and generally 
falls within the experimental error for the data. The wall pressure and density are in agree- 
ment within the experimental error. On the basis of experimental error, the pitot pressure 
shows somewhat more variation than is expected; this variation may be due to an increase in 
the flow velocity behind the interface, as reported in reference 4. The rapid increase in 
density and decrease in temperature early in the flow is, in part, a result of the response 
time (2.5 ps) of the system but may also be caused by mixing at the interface. It should be 
recalled that the measured temperature is a vibrational temperature and the temperature inferred 
from the density and wall pressure is a translational temperature; the difference between them 
is probably an indication of vibrational nonequilibrium. 

A second example, run 3679, of the time history of the flow is shown in figure 9. The 
initial conditions upstream of the secondary diaphragm for this run are about the same as for 
the previous example; however, no helium was added in the acceleration section. Thus the 
flow has been expanded to a very low pressure, and there is no detectable region of uniform 
flow conditions. The pressure and density measurements are nearly ramp functions with time, 
and the vibrational temperature takes about 100 jus to come to equilibrium with the static 
temperature instead of the 25 ps in the previous example. This longer relaxation time is 
qualitatively what would be expected if the flow were not in vibrational equilibrium during 
the early part of the flow. 


Data Trends 

As a means of obtaining data trends, tests have been grouped according to test gas and 
shock speed, and the results are plotted as a function of Vj. Because the pressure at which 
the first diaphragm ruptured was not the same for each run, there is a small scatter in Vg j 
(about 1.5 percent). This effect was eliminated by plotting the conditions for an average 
V§ j. In order to correct for the differences introduced between the measured and theoretical 
predictions by this averaging process, the measured data were scaled using the RSEE and RSFE 
models. Thus the location of the data relative to the RSEE and RSFE predictions is the same 
for both the true and the averaged test condition. This procedure is described in detail in 
reference 1. 

The variations of p, T, p w , and p t with interface speed for a series of tests with 
air as a test gas are shown in figure 10. The circular and square symbols denote the averages 
of the measured data over a period of 50 ps in t^ from t<. = 50 ps to t^ = 100 ps. 

The triangular symbol denotes the corresponding value inferred from equation (8) or (9). In 
addition, if the measured value can be determined over the first 50 ps of the flow and if it 
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is significantly different from the average over the second 50 ps, the range of these two 
averages is indicated by a bar. The predicted curves were obtained by using pj, Tj, Vg 
Vj and the computer code described in reference 17. 

The measured densities, figure 10(a), always lie above the RSFE model and, with but 
two exceptions, are below the predicted value of any of the other models. The values of 
density inferred from equation (8), p w> and T V; Q 2 lie near or below the measured value 
and show better agreement with the measured values for the higher velocity cases. The tests 
made with the diaphragm opener show much better agreement between the inferred and mea- 
sured values although this could be a velocity effect rather than a result of the diaphragm 
opener. 

The variation of TyC^ i s shown in figure 10(b). The theoretical predictions here refer 
to the translational temperature, not to the vibrational temperature. For this case the vibra- 
tional temperatures for the ISFE and RSFE models are 3169 and 4985 K, respectively. With 
but two exceptions the inferred temperature (obtained from eq. (8), p w , and p) lies well 
below the measured value. This indicates that the flow is not in vibrational equilibrium. 

The measured wall pressures, figure 10(c), generally fall between the ISFE and ISEE 
predictions. The values inferred from equation (8) are generally higher, with the best agree- 
ment between the measured and inferred values again occurring for the higher interface speeds 
and the diaphragm opener runs. 

The measured pitot pressures, figure 10(d), show almost no systematic variation with 
interface speed. The values inferred from equation (9) are generally above the measured 
values, with the only two exceptions being two diaphragm opener tests. In general, the data 
lie between the frozen predictions. 

The data trend as a function of interface speed for a series of tests with as the 
test gas is shown in figure 11. Because the accuracy of the density and temperature mea- 
surement depends on the amount of in the light beam, these tests were made at a lower 
density than were the air runs. 

The density measurements, figure 11(a), lie about the ISFE prediction and show good 

agreement between the measured and inferred values. Also shown in figure 11(a) are two 

runs at Vj of 6250 and 6350 m/s which were made without adding any helium to the 
acceleration chamber and thus show large changes in the measurements since there is no 
period of steady flow. 

The vibrational temperature, figure 1 1 (b), shows some good agreement between the mea- 
sured and inferred values for the second 50 ps average but also shows strong relaxation with 

time for the first 50 ps. This indicates that while the vibrational temperature is far from the 
translational temperature early in the flow, it is approaching equilibrium with time. 
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The wall pressure, figure 11(e), falls just below the ISPG model and again shows good 
agreement between the measured and inferred values. 

The pitot pressure, figure 11(d), lies about the ISFE prediction and shows fair agreement 
between the measured and inferred values for most cases. 


These results for air and indicate that none of the theoretical models examined pro- 
vides a means of predicting the state of expansion tube test flow over the range of conditions 
tested. One possible explanation for this is the fact that the flow is highly sensitive to small 
levels of dissociation. This sensitivity is illustrated by figure 11(d). The mass fraction at 
the test section for the ISFE model is 0.94; thus, if the flow were described by the model, 
the pitot pressure inferred from equation (9) should be 6 percent less than the measured pitot 
pressure. However, the predicted values obtained from the ISEE and ISFE models differ by' 
more than 250 percent. 


In order to accentuate the effects of dissociation, a series of tests were made with 

at a reduced pressure as the test gas. This resulted in a stronger primary shock and higher 

degree of dissociation upstream of the expansion fan. The ratio A for these tests is shown 

in figure 12, where it has been plotted as a function of the stagnation enthalpy (calculated 

from ref. 17) behind the reflected shock at the secondary diaphragm. The mass fraction F 

has been assumed to be unity, and equation (9) has been used to calculate p ; thus A 

t,c 

becomes 


A = 


P °2 

p o 2 + p o 


(14) 


Consequently, A should be unity if the flow is in thermodynamic equilibrium and less than 
unity for any degree of nonequilibrium. Figure 12 shows A to be in good agreement with 
the RSFE model for H 6 greater than 10 7 m 2 /s 2 ; however for H g less than this, the 
differences between the models are less than the experimental errors in A. The reader should 
not assume that figure 12 shows the RSFE model to be the proper model for these flow 
conditions; it only shows that to a good approximation the degree of dissociation in the test 
flow is predicted by the stagnation conditions upstream of the secondary diaphragm. 

A similar plot for a group of air tests is shown in figure 13. Here the data have been 
plotted as a function of H^ because the dr + a trend is not in agreement with the RSFE 
model. The data fall into two distinct groups — those with the diaphragm opener and those 
without. When the diaphragm opener was used, the ratio was less than unity and was in good 
agreement with the ISFE model; however, when the diaphragm opener was not used, the ratio 
was greater than unity. This result would imply that the mass fraction of in the test 
flow was greater than that in the initial test-gas mixture. Similar results were reported in 
reference 1, and the authors still have no explanation for this effect. 
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CONCLUDING REMARKS 


Simultaneous time-resolved measurements of temperature, density, pitot pressure, and 
wall pressure in both air and O2 test gases have been obtained in the Langley pilot model 
expansion tube. These tests show that nonequilibrium chemical effects are present in the 
test flow and significantly affect the calculation of the flow condition. These effects must 
be accounted for from the primary shock since calculations have shown the residence time 
behind the primary shock is too short to allow the flow to reach equilibrium before entering 
the expansion fan; measurements have shown that the degree of dissociation in the test gas 
is determined by conditions upstream of the second diaphragm for some test conditions. 

The measured vibrational temperatures were higher than both the equilibrium predictions 
and the translational temperatures inferred from the measured density and wall pressure. The 
vibrational temperature showed significant relaxation during the test flow, thereby demon- 
strating that vibrational nonequilibrium effects were most significant early in the flow. 

The use of splitter plates to reduce the boundary-layer effects on the density and tem- 
perature measurement showed that the discrepancy between the measured pitot pressure and 
the value inferred from the measured density and interface velocity was not a result of the 
measurements being made across the tube boundary layer. 

The use of an electromagnetic device to open the secondary diaphragm prior to the 
arrival of the primary shock resulted in an improvement in the agreement between the mea- 
sured and inferred pitot pressure. Because of the small numbers of tests with the opener, this 
result is not conclusive; however, it does indicate that the effects of the opening of the secon- 
dary diaphragm may have a strong influence on the test flow. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
October 22, 1975 
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TABLE 1.- TEST CONDITIONS 


Run 

Pp Pa 

PlO* Pa 

V S, 1’ m/s 

v j> m/s 

Test gas 



Diaphragm opener 



3646 

2933 

26.7 

2718 

5746 

Air 

3647 

2933 

13.3 

2942 

5499 

Air 

3648 

2933 

6.67 

2873 

5390 

Air 

3649 

2933 

2.67 

2808 

5806 

Air 



Splitter plates 



3660 

2933 

o 

Tj- 

— 

2855 

5034 

Air 

3661 

2933 

133 

2929 

5240 

Air 

3662 

2933 

133 

2917 

4938 

Air 

3663 

2933 

133 

2779 

4961 

Air 

3664 

2933 

133 

2996 

5053 

Air 

3665 

2933 

66.7 

2872 

5667 

Air 

3666 

2933 

66.7 

2886 

5291 

Air 

3667 

2933 

40 

2878 

5451 

Air 

3668 

2933 

40 

2874 

5046 

Air 

3669 

2933 

66.7 

2911 

5083 

Air 

3670 

2933 

66.7 

2860 

4980 

Air 

3674 

2933 

26.7 

2961 

5678 

°2 

3675 

2933 

13.3 

2912 

5730 

°2 

3676 

2933 

13.3 

2955 

6101 

°2 

3677 

2933 

13.3 

2806 

6014 

02 

3678 

2933 

6.67 

2832 

6112 

°2 

3679 

2933 

(a) 

2865 

6247 

o 2 

3680 

2933 

(a) 

2858 

6354 

o 2 

3683 

1467 

13.3 

3310 

6458 

°2 

3684 

667 

13.3 

3607 

6236 

o 2 

3685 

333 

13.3 

4002 

6212 

.0 2 

3686 

1467 

13.3 

3331 

6361 

o 2 

3687 

667 

13.3 

3767 

6347 

°2 

3688 

333 

13.3 

3834 

6325 

°2 

3690 

2933 

6.67 

2957 

6003 

°2 


a Less than 1 Pa. 




SECTION SECTION SECTION 

(a) Expansion tube flow cycle showing pertinent flow regions. 



(b) Expansion flow cycle at region of second diaphragm, with effects 

of diaphragm considered. 

Figure 1.- Distance-time diagrams. 








VACUUM SYSTEM 



(a) Test section without splitter plates. 


VACUUM SYSTEM 



(b) Test section with splitter plates. 

Figure 2.- Sketch of test section and lamp assembly. 
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Figure 4.- Installation of coil in expansion tube. 



(g) t = 140 jus. 


(h) t = 190 fis. 


(i) t = 240 [xs. 

L-75-237 

Figure 5.- Diaphragm opening process for intervals of time t after trigger. 





Figure 6.- Temperature dependence of absorption coefficients 
and their ratio (from ref. 14). 



Figure 7.- Mole fraction of O upon entrance into expansion fan. . 
Pj = 2932 Pa; Tj = 300 K; V g ^ = 2956 m/s; test gas: Og. 
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Figure 8.- Measured and inferred flow quantities for 
run 3676 of table 1. pj = 2933 Pa; Tj = 300 K; 
V S,1 = 2955 m / s ; v i = 6101 m /s; test gas: 0 2 . 



Figure 9.- Measured and inferred flow quantities for 
run 3679 of table 1. pj = 2933 Pa; Tj = 300 K; 
V Sj i - 2865 m / S ; V i = 6247 m/s; test gas: O 2> 



Experiment 

O Measured with splitter plates 
E2 Measured with diaphragm opener 



4.6 4.8 5.0 5.2 5.4 5.6 5.8 6.0 

V j , km/s 


(a) Density. 

Figure 10.- Comparison of measured, inferred, and predicted test data 
for air test gas. V g j = 2880 m/s; p 1 = 2932 Pa; Tj = 300 K. 
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Experiment 

O Measured with splitter plates 
Eg Measured with diaphragm opener 
A Inferred from measurements 
and eq. (8) and (9) 

Model predictions (ref. 17) 

RSEE 

RSFE 

ISEE 

ISFE 




5.2 5.4 5.6 5.8 6.C 

Vj , km/s 


(c) Wall pressure. 
Figure 10.- Continued. 



Experiment 

O Measured with splitter plates 
^ Measured with diaphragm opener 
A Inferred from measurements 
and eq. (8) and (9) 


Model predictions (ref. 17) 
— RSEE 



(d) Pitot pressure. 
Figure 10.- Concluded. 
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Experiment 

O Measured with splitter plates 
A Inferred from measurements 
and eq. (8) and (9) 

Model predictions (ref. 17) 
RSEE 



6.0 6.2 6.4 

V., km/s 


(b) Temperature. 
Figure 11.- Continued. 



Experiment 

O Measured with splitter plates 
A Inferred from measurements 
and eq. (8) and (9) 



(c) Wall pressure. 
Figure 11.- Continued. 
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Experiment 

Measured with splitter plates 
Inferred from measurements 
and eq. (8) and (9) 


2 


Model predictions (ref. 17) 

RSEE 

ISEE 

ISFE 

ISPG 



(d) Pitot pressure. 
Figure 11.- Concluded. 







O Splitter plate 
□ Diaphragm opener 



Figure 13.- Variation of A with H 2 for air test gas. 
Pj = 2933 Pa; Tj = 300 K. 
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